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We describe a heat capacity measurement system for the study of the field-orientation dependence 
for temperatures down to 50 mK. A "vector magnet" combined with a mechanical rotator for the 
dewar enables the rotation of the magnetic field without mechanical heating in the cryostat by 
friction. High reproducibility of the field direction, as well as an angular resolution of better than 
0.01°, is obtained. This system is applicable to other kinds of measurements which require a large 
sample space or an adiabatic sample environment, and can also be used with multiple refrigerator 
inserts interchangeably. 



I. INTRODUCTION 



The most important advances in superconductivity re- 
search over the past two decades have been the discover- 
ies and studies of superconductors in which strong elec- 
tronic correlations play crucial roles in both the nor- 
mal and superconducting (SC) states. Thus, the so- 
called "unconventional" superconductivity, realized in 
high-Tc cuprates, organics, heavy-fermion intermetallic 
compounds, ruthenate^ etc., has become one of the most 
actively studied topics of modern-day condensed mat- 
ter physics. Since the anisotropic structure of the SC 
gap is directly connected with the SC order parameter 
and is closely related to the anisotropy of the pairing in- 
teraction, experimental determination of the gap struc- 
ture is essential for the identification of the mechanism of 
superconductivity, especially among the unconventional 
superconductors. The field-orientation dependent spe- 
cific heat, which is a direct measure of the quasiparti- 
cle density of states, is a powerful probe of the SC gap 
structureiSi^ 

Such an experiment requires: (1) reaching a very 
low temperature compared with the SC transition tem- 
perature Tc, (2) high angular resolution of exact field 
alignment with respect to crystallographic axes, and (3) 
high- resolution calorimetry in magnetic fields. To sat- 
isfy condition (1), we utilized a commercial dilution re- 
frigerator. Unconventional superconductors often have 
low-dimensional electronic structures and exhibit strong 
anisotropy in the upper critical field i?c2 (Ref. 0). Thus 
the angular resolution of better than 0.1° is often re- 
quired for the measurements. To control the field direc- 
tion, the rotation of a sample with a geared goniome- 
ter has been widely used. A single-axis goniometer, ^'^ a 
low friction single-axis goniometer^ and a double-axis 



goniometer" 



8.9.10 



have been used. As a single-axis go- 



niometer with a very compact and simple structure, a 
piezoelectrically driven goniometer is also in use>ii For 
rotation in liquid-helium mixtures, a miniature rotatable 



vacuum cell has been developed. However, many of the 
geared goniometers cannot avoid the heat generated by 
mechanical friction and suffer from poor angular repro- 
ducibility caused by the backlash of the gears. To achieve 
the condition (2) and solve these problems, we developed 
a new system in which we electrically rotate the magnetic 
field with a "vector magnet" and mechanically rotate this 
SC magnet with a one-circle goniometer against a fixed 
dilution refrigerator. 

The specific heat measurements in magnetic fields (3) 
have been carried out by a relaxation method^"^ with an 
automated small sample calorimeterii^ii^ from 50 mK to 
8 K using a dilution refrigerator >iSiii In the measure- 
ments, we often investigate the field amplitude \H\ and 
field orientation (/>, 9 dependence of the specific heat at 
constant temperatures and need to know the absolute 
value of the specific heat. In this article, we describe the 
construction, performance, and advantage of this system. 



II. DESIGN AND CONSTRUCTION 

For controlling the field orientation at low tempera- 
ture, we built a system with two orthogonally arranged 
SC magnets (Cryomagnetics, Inc., Oak Ridge TN, model 
VSC-3050) which generates the field H with a horizon- 
tal field Hr and a vertical field Hz to control the polar 
angle 9 as shown in the inset of Fig. ^ A split pair 
and a solenoidal SC magnet provide of up to 5 T 
and Hz of up to 3 T, respectively, with the homogene- 
ity of ± 0.1% over 10 mm diameter spherical volume. 
This vector magnet allows control of the field orientation 
9 = aicta.n{Hz/Hr) with a precision A9 = 0.001°. 

The two SC magnets are installed at the bottom of 
a dewar seating on a goniometer (HUBER Diffraktion- 
stecknik GmbH, Rimsting Germany, one-circle goniome- 
ter model 420) which controls the azimuthal angle of 
the field direction. The assembly can be rotated relative 
to a dilution refrigerator (Oxford Instruments, Inc., Wit- 
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FIG. 1: Schematic diagram of the specific heat measurement 
system with a vector magnet. 



ney, Oxon, UK Kclvinox 25) which is fixed, with a mini- 
mum step of Acf) — 0.001°. The goniometer can rotate the 
two SC magnets and the dewar with the magnetic field 
applied either in a persistent-current mode or in a cur- 
rent feeding mode. We used long enough magnet cables 
to allow rotation of ±360°. Controlling the two magnets 
and the goniometer, we can rotate and align the field H 
continuously to any direction. For example, for the study 
of a system with a two-dimensional (2D) electronic struc- 
ture, we can rotate the field continuously within the 2D 
plane with a misalignment no greater than 0.01° from the 
plane. Advantages of the field rotating system with the 
vector magnet and the goniometer, compared with a sam- 
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applicability to different refrigerator inserts interchange- 
ably, (3) no mechanical heating by friction at low temper- 
atures except for the eddy current heating by variations 
of the magnetic field, and (4) good reproducibility of the 
field orientation. 

Details of the addenda of the calorimeter are schemat- 
ically shown in Fig. [3 The thermometer is made of a 1.8 
kfl Ru02 thick-film resistor (KOA, RK73KlEJ).i« To 
maximize the thermal response, as well as to minimize 
the heat capacity of addenda, the alumina substrate of 
the Ru02 thick-film resistor is polished and thinned down 
to the size of 1.0 x 0.5 x 0.1 mm^. The solder on the ter- 
minal of the Ru02 thick-film resistor is also shaved to 
eliminate a change of addenda heat capacity at the SC 
transition of the solder. The heater is made of a 350 
strain gauge (Kyowa, KFG-2-350-Cl-16N2C2)J.a which is 
cut out to the size 2.5 x 2.5 x 0.05 mm"^. The thermome- 
ter and heater are glued directly on one side of a thin 
foil of oxygen-free high-conductivity copper with 5N pu- 
rity, 3.0 x 3.0 X 0.025 mm'^, using a minimal amount of 
epoxy (Emerson and Cuming, Stycast 1266). The sam- 
ple is placed on the other side of the copper foil with 
a small amount of grease (Apiezon Products, N grease). 
The electrical leads are monofilament niobium-titanium 
wires, 67 /im in diameter, with a SC transition temper- 
ature Tc = 9.3 K and upper critical field /ioi?c2 = 12 
T. We use this wire in the SC state in order to elimi- 
nate the parasitic heating as well as to suppress the heat 
leak. These wires are attached to the thermometer and 
heater by heat-cure-type silver epoxy (Dupont, 6838). 
As shown in Fig. 13 a platinum wire with the diameter of 
10 /im is connected with normal-cure- type silver epoxy 
(Dupont, 4922N) as a thermal link between the addenda 
and the copper base, which serves as a thermal bath. The 
calorimeter is covered with a copper foil as a radiation 
shield. 

Heat capacity measurements are carried out with a re- 
laxation method^'^ as follows. When we apply a constant 
power Pq to the addenda, the temperature increases from 
the base temperature Tq to Tq + AT at a steady state. 
The thermal conductance K of the wires is evaluated as 
K = Pq/AT. By switching off the power, the tempera- 
ture relaxes from Tq + AT down to Tq according to 



T{t) ^To + ATexp {-t/n). 



(1) 



FIG. 2: Schematic drawing of the calorimeter with an ad- 
denda and a copper base for the relaxation method. 



We determine the relaxation time ri from the fitting to 
a relaxation curve in the range of < t < lOri. We typi- 
cally obtain ri ~ 25 s in the heat capacity measurement 
of aluminum and Sr2Ru04 as described below. Thus, 
the total heat capacity can be calculated from C = Kti . 
The simple model given above yields good results in most 
cases, but if the thermal conductance between the sample 
and addenda cannot be made sufficiently large compared 
with that of the thermal link, the relaxation curve is char- 
acterized by a large initial slope compared to the rest of 
the decay, which is called the T2 effect i^SiSi In that case, 
the relaxation curve is better described by 



T{t) ^Tq + A exp(-t/Ti) + B exp(-t/T2) 



(2) 
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FIG. 3: Magnetoresistance AR/R = [R{H,T) - 
R{0, T)]/R{0, T) vs magnetic field for a Ru02 thick film resis- 
tor at several temperatures. Inset: Temperature dependence 
of resistance of a Ru02 thick film resistor in zero field. 
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FIG. 4: Magnetoresistance AR/R = [R{H,T) - 
R{0,T)]/ R{0,T) vs magnetic field orientation for a Ru02 
thick-film resistor at T = 0.14 K and fioH = 2 T. The po- 
lar angle 6 is defined as the angle between the field direction 
and the axis perpendicular to the Ru02 film plane, while the 
azimuthal angle is defined as the angle from the current 
direction. 



The total heat capacity is described as C = K{Ati + 
Bti)I AT with A + B = AT. Improving the thermal 
contact between the assembly of addenda and a sam- 
ple results in very small effect: We typically obtain 
T2/T1 < 0.05 and BjA < 0.05 in the range of 50 mK 
< T < 8 K. Thus the total heat capacity is approximated 
as C ~ Kti-^. Since precise measurements of short T2 
are difficult with nonlinear fits of relaxation curves and 
« 1 is obtained in our case, we use this approximated 
equation. Therefore, the accuracy of heat capacity mea- 
surement with a relaxation method is well assured in this 
apparatus. 

For a relaxation curve measurement, which requires 
fast readings of the thermometer resistance, we use 
a phase-sensitive detector (Stanford Research, model 
SR830) with the measurement frequency of 439 Hz and 
a 30 ms time constant. To give a heat pulse and measure 
the heater power, a dc current source (Yokogawa, model 
7651) and a nano-voltmeter (Keithley, model 2182) are 
connected to the addenda heater. In order to accurately 
measure the specific heat, it is very important to main- 
tain a constant temperature of the thermal bath. For 
this purpose, the temperature of the copper base is mea- 
sured with a thermometer (a 1.8 kO Ru02 thick film 
resistor) and an ac resistance bridge (RV-Elektroniikka, 
model AVS-47) using the excitation frequency of 15 Hz, 
while it is controlled with a heater (strain gauges with 
350 Vl) and a temperature controller (RV-Elektroniikka, 
model TS-530). To avoid noise pick up and rf heating 
of the sensors, the signal wires are filtered at the top 
of the cryostat at room temperature by means of sim- 



ple low-pass LC (inductance-capacitance) filters with a 
cut off frequency of 10 kHz. Consequently, we achieved 
the thermometry with low noise and high stability at low 
temperatures. 



III. CALIBRATION 

The calorimetry in magnetic field requires appropri- 
ate calibrations of the resistance thermometer for its 
magnetoresistance z4i?/i? = [R{H,T)-R{Q,T)]/ R{Q,T). 
In addition, the examination of anisotropy of its mag- 
netoresistance is important in the study of the field- 
orientation dependence of the specific heat. The RUO2 
thick- film resistor with a nominal resistance of 1.8 kfi 
(KOA, RK73K1EJ) was adopted as the resistance ther- 
mometer of the calorimetry in magnetic field because 
of its advantages of high sensitivity at low temperature, 
small heat capacity, fast thermal response time, and low 
costi^^i^3,24,25,26,27,28,29 the cahbration of the ther- 

mometer in zero field, commercially calibrated RUO2 
thermometer (Scientific Instruments, model RO600A2) 
was used in the range from 50 mK to 8 K. The inset in 
Fig. 13 shows the temperature dependence of resistance 
of the addenda thermometer in zero field. It is well 
represented by a characteristic temperature dependence 
of the electrical conduction mechanism due to variable 
range hopping in three-dimensional Anderson localized 
states^": R = Ro exp(To/T)i/4. 

Field-strength and field-orientation dependences of the 
Ru02 thick film resistor were measured at a constant 
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TABLE I: Comparison of normal-state and SC-state param- 
eters of aluminum with previous report. 
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The coefficients am,n are determined from the data taken 
from the field sweeps at 12 temperatures. A part of the 
results of fitting made in the range of 50 mK < T < 8 K, 
T < noH < 4.25 T are shown in Fig.El 



FIG. 5: Temperature dependence of the specific heat divided 
by temperature C/T of aluminum in several magnetic fields. 
Inset: Field dependence of C/T at T = 0.10 K. 



temperature. To maintain a constant temperature, each 
measurement was carried out in the persistent-current 
mode of the SC magnet after very slow field sweep- 
ing and establishing a thermal equilibrium state. Fig- 
ure 131 shows the field-strength dependence of the mag- 
netoresistance at several temperatures. At low temper- 
atures, the magnetoresistance is negative at low fields 
but becomes positive at higher fields, consistent with 
previous reportsi2^iSS*2L2& Such behavior is attributable 
to disorder effects (weak localization) at low tempera- 
tures, which produce negative magnetoresistance in the 
low-field region. '^^•'^^ At high fields, the field dependence 
is nearly proportional to the square root of H, which 
originates from weak localization with electron-electron 
scattering iSi 

Figure0|shows the anisotropy of the magnetoresistance 
at 0.14 K, investigated in three ways: The azimuthal an- 
gle (j) dependence within the plane of the Ru02 film, the 
polar angle 9 dependence from the current direction, and 
the polar angle 6 dependence perpendicular to the cur- 
rent direction. As is clear from Fig.^ no significant field- 
orientation dependence against the Ru02 chip and cur- 
rent directions was found. This feature provides an ex- 
cellent advantage to the thermometry in magnetic fields. 
The information in Fig. ^ explicitly represents full sets 
of data on the field-orientation dependence of a Ru02 
for the first time, which is consistent with comments in 



previous reports 



22.25 



Since the dependence of 9 and 



can be neglected, the resistance R{T, H) is expressed in 
terms of a 2D-polynomial function of temperature T and 



IV. PERFORMANCE 

In order to demonstrate its performance, we examine 
the accuracy of the measured specific heat of aluminum 
(Al), a conventional Type-I superconductor with Tc = 
1.17 K. We also inspect the field-orientation dependence 
of the heat capacity of the addenda. A polycrystalline 
aluminum sample with a purity of 6N was cut from the 
polycrystalline rod, to a size of 3.9 x 3.9 x 1.9 mm'^ with 
77.9 mg. The specific heat C of the sample was evaluated 
from the total heat capacity Ctotai after subtraction of 
the addenda heat capacity Caddenda. The heat capacity 
of the addenda was separately measured and is described 
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FIG. 6: Field-orientation dependence of the total specific 
heat Ctotai = C + Caddenda divided by temperature of the 
aluminum in the normal state at T = 0.10 K and fioH — 20 
mT. 



5 



as Caddonda = (« + a'H^)/T^ +jT + (3T^ with a = 
1.03 X 10"'' nJK, a' = 4.68 x lO^^ ^JK/T^, 7 = 1.15 x 
10-1 mJ/K2, and P = 4.36 x 10"^ ^J/K". 

Figure shows the temperature dependence of C/T 
of the aluminum sample in magnetic fields. At fJ-oH — 
T and 7.5 mT, the SC transitions were observed at 
Tc = 1.16 K and 0.62 K, respectively. At ^iqH = 20 
mT, the normal state persists down to 50 mK. From the 
specific heat in the normal state, we can estimate the 
electronic specific heat coefficient 7 = 1.37 mJ/K^mol 
and the Debye temperature Or, — 420 K from the fitting 
to C = 7T + ^T3,/3 = {12/ 5)7:'^ R6^^ below 4 K. In zero 
field, the specific heat jump AC/jTc — 1.37 was observed 
with a SC transition at Tc — 1.16 K. These results are 
in excellent agreement with a previous report by an adi- 
abatic method using a sample weighting 212.64 gj^ as 
summarized in Tabled The temperature dependence of 
the specific heat is fully consistent with theoretical cal- 
culation based on the Barden-Cooper-Schrieffer theorjiSi 
in all temperature range. This examination assures the 
high accuracy and high resolution in the specific heat 
measurement with this system. The requirement of the 
entropy balance suggests that the transition under the 
magnetic field {iioH = 7.5 mT) is of the first order with 
a latent heat of L = T^AS = 0.091 mJ/mol. 

The inset of Fig. [S] shows the field dependence of C /T 
at T = 0.10 K. C/T is nearly equal to the electronic con- 
tribution Cc/T oc density of state (DOS) at low temper- 
atures. Since a fully gapped SC gap opens and removes 
the DOS at the Fermi surface, C/T remains nearly zero 
at low fields. With an increasing magnetic field, the DOS 
increases because of the normal state region in the inter- 
mediate state, and Ce/T oc DOS recovers up to 7 at the 
critical field He = 10.5 mT. This value is also consistent 
with the previous report^ and guarantees the accuracy 
of the applied magnetic field amplitude generating by the 
vector magnet. 

Figure shows the field-orientation dependence of the 
total specific heat divided by temperature Ctotai/^ with 
the aluminum sample in the normal state at T = 0.10 
K and hqH = 20 mT. The fraction of the addenda heat 
capacity is Caddcnda/C'totai = 0.08. The definitions of 
9 and (j) are with respect to the addenda thermometer 
as shown in Fig. 0] Since the polycrystalline aluminum 
in the normal state does not give the field-orientation 
dependence of the total heat capacity, these results re- 
veal that the addenda heat capacity has totally negligible 
field-orientation dependence. 



V. APPLICATION EXAMPLE 

As an application example, we describe an exper- 
iment for determination of the SC gap structure of 
Sr2Ru04 from the field-orientation dependence of the 
specific heat. The layered perovskite Sr2Ru04 is 
an unconventional Type-II superconductor with strong 
quasi-two dimensionalityj^ Its superconductivity ex- 
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FIG. 7: (a) The polar angle 9' dependence of the electronic 
specific heat divided by temperature Ce/T of Sr2Ru04. (b) 
The azimuthal angle (p' dependence of Ce/T. (c) Field de- 
pendence of Cc/T in the vicinity of the upper critical field 

Hc2. 



hibits pronounced unconventional features attributable 
to spin-triplet p-wave pairing with time-reversal symme- 
try breaking, or "chirality" , associated with the Cooper- 
pair orbital moment. ^ In spite of a large number of ex- 
periments and theories on its SC gap structure and spin- 
triplet superconductivity, the details of the gap structure 
still remain to be clarified. 

Single crystals of Sr2Ru04 were grown by a floating- 
zone method in an infrared image furnace»2& The sample 
with Tc = 1.48 K, close to the estimated value for impu- 
rity and defect-free specimen (Tco — 1.50 K), was chosen 
for detailed study. This crystal was cut and cleaved from 
a single crystalline rod, to a size of 2.8 x 4.8 mm^ in the ah 
plane and 0.50 mm along the c axis, and weighs 39.8 mg. 
An x-ray rocking curve of the sample shows the charac- 
teristics of a single crystal of high quality; the diffraction 
peak width [full width at half maximum (FWHM)] be- 
ing comparable to that for a Si crystal (with a FWHM 
of 0.06°) limited by the instrumental resolution of our 
diffractometer. The directions of the tetragonal crystal- 
lographic axes of the sample were determined by x-ray 
Laue pictures. The side of the crystal was intentionally 
misaligned from the [110] axis by 16° to distinguish the 
essential anisotropy related to the unconventional super- 
conductivity from the extrinsic anisotropy due to surface 
effects. 

In this experiment, capability of rotating H within the 
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FIG. 8: (a) The in-plane field-orientation dependence of 
the electronic specific heat divided by temperature Ce/T of 
Sr2Ru04 at ^oH = 0.6 T and several temperatures, (b) The 
phase diagram of Sr2Ru04 for H \\ [100], based on specific 
heat measurements. Hc2 is the upper critical field and crosses 
indicate where data presented in (a) are taken. 



Ru02 plane with high accuracy is very important be- 
cause of strong quasi-2D electronic structure with cylin- 
drical Fermi surfaces. A slight field misalignment causes 
extrinsic two-fold anisotropy of the field-orientation de- 
pendent specific heat due to a large Hc2 anisotropy 
{Hc2\\ab/Hc2\\c ~ 20).'^''-^^ For specifying the direction of 
the applied magnetic field with respect to the crystallo- 
graphic axes, we shall introduce the polar angle 6', for 
which 9' = 0° corresponds to the [001] direction, and 
the azimuthal angle 0', for which cj)' = 0° with 6' = 90° 
corresponds to the [100] direction. It is noted that the 
tetragonal symmetry of the crystal structure is conserved 
down to temperatures as low as 110 mKi^ 

The field alignment was carried out with the polar an- 
gle 9' dependence of specific heat at several azimuthal 
angle (p' . With the choice of an appropriate field ampli- 
tude Hc2\\ab ^ H Hc2\\c, thc SC statc emerges only for 
the magnetic field nearly parallel to the Ru02 plane, and 
Cc/T cx DOS decreases in the SC state because DOS is an 
increasing function of H/Hc2 at low temperature, just as 
shown in the inset of Fig. [S] for aluminum. Since Hc2{9') 
takes a maximum at H \\ Ru02 plane, Cc/T takes a min- 
imum. Figure [7{a) shows the 9' dependence of Cc/T of 
Sr2Ru04 at T = 0.11 K and /io-ff = 1.38 T. Note that 
the result clearly resolves that IjI.oHc2{9') > 1.38 T only 
within A9' ~ ±1.3° with high angular precision. From 
this procedure at several azimuthal angles, we can ob- 
tain information to align H along the Ru02 plane with 



an accuracy better than 0.05°. Since the 9' dependence of 
Cc/T is symmetric against 9' = 90°, the effect of vortex 
pinning in the 9' sweep was negligibly small for Sr2Ru04. 

The in-plane field alignment relative to the crystal- 
lographic axis 4>' was carried out, based on the known 
in-plane anisotropy of Hc2- i^c2||[iio] > -f^c2||[ioo] at low 
temperaturesi^Si^ Figure (Tfb) shows the 0' dependence 
of Cc/T at T = 0.51 K and hqH = 1.39 T. A sinusoidal 
four-fold angular variation of Cc{4>)/T is observed. Fig- 
ure[7Jc) shows the field dependence of Cc/T at H \\ [110] 
and H II [100] at T = 0.51 K. In the vicinity of i7c2, Cc/T 
takes maxima in magnetic fields parallel to the [110] di- 
rections because ^fc2||[iio] > ^c2||[ioo]- From these re- 
sults, we can assign the directions of the crystallographic 
axes in the Ru02 plane with an accuracy better than 
2.0°. 

In the mixed state, the quasiparticle (QP) energy spec- 
trum Ek is affected and becomes E'j^ = — Suj by 
thc Doppler shift 6lu ~ hk ■ Vs, where Vg is the super- 
fluid velocity around the vortices and hk is thc QP mo- 
mentum. In the case of 5lu > A{k), this energy shift 
gives rise to a finite DOS at the Fermi level3 Since the 
superfluid velocity Vg is perpendicular to the magnetic 
field H, the Doppler shift Scu becomes zero for k \\ H. 
Thus, the generation of QPs is suppressed for H \\ gap 
minima directions and Cc/T oc DOS takes minima for 
these directionsi^ Consequently, the field-orientation de- 
pendence of the specific heat directly reveal the directions 
of the gap minima. 

Figure ISJa) shows the field-orientation dependence of 
the specific heat of Sr2Ru04. The absence of a two- 
fold oscillatory component in the raw data guarantees 
that the in-plane field alignment is accurate during the 
azimuthal- angle (/>' rotation. At iiqH — 0.6 T and low 
temperatures, shown in Fig. |HIb) as crosses, a non- 
sinusoidal four-fold angular variation approximated as 
6Cc{(t))/T oc 2|sin20'| - 1 is observed. Since 5Cc{(l))/T 
due to the in-plane Hc2 anisotropy would have an in- 
verse phase below 1.35 T and the amplitude of the an- 
gular variation in Fig. IHl^a) becomes smaller with an in- 
crease of temperature toward Tc, we conclude that the 
4-fold oscillations at a low temperature originate from 
the anisotropy in the SC gap. 

We have deduced the SC gap structure in the inset 
of Fig. ISl^b) with a gap minimum along the [100] direc- 
tion, because Cc{(t>)/T takes a minimum to the [100] di- 
rection and shows in-plane four-fold symmetry. The SC 
gap function is well described as A{k) = Ao{sm^akx + 
sin^afcy)^/^ with the gap minima; this SC gap structure 
is promising for the active SC band of Sr2Ru04. The 
details of this study are published elsewhere^i 
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